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The magnetoelectric behavior of BiFeOa has been explored on the basis of accurate density func- 
tional calculations. We are able to predict structural, electronic, magnetic, and ferroelectric prop- 
erties of BiFeOs correctly without including any strong correlation effect in the calculation. Unlike 
earlier calculations, the equilibrium structural parameters are found to be in very good agree- 
ment with the experimental findings. In particular, the present calculation correctly reproduced 
experimentally-observed elongation of cubic perovskite-like lattice along the [111] direction. At 
high pressure we predicted a pressure-induced structural transition from rhombohedral (i?3c) to an 
orthorhombic (Prima) structure. The total energy calculations at expanded lattice show two lower 
energy ferroelectric phases (with monoclinic Cm and tetragonal Pimm structures), closer in energy 
to the ground state phase. Spin-polarized band-structure calculations show that BiFeOs will be 
an insulator in A- and G-type antiferromagnetic phases and a metal in C-type antiferromagnetic, 
ferromagnetic configurations, and in the nonmagnetic state. Chemical bonding in BiFeOs has been 
analyzed using partial density-of-states, charge density, charge transfer, electron localization func- 
tion. Born-effective-charge tensor, and crystal orbital Hamiltonian population analyses. Our electron 
localization function analysis shows that stereochemically active lone-pair electrons are present at 
the Bi sites which are responsible for displacements of the Bi atoms from the centro-symmetric to 
the noncentrosymmetric structure and hence the ferroelectricity. A large ferroelectric polarization 
of 88.7/iC/cm^ is predicted in accordance with recent experimental findings, but differing by an 
order of magnitude from earlier experimental values. The strong spontaneous polarization is related 
to the large values of the Born-effective-charges at the Bi sites along with their large displacement 
along the [111] direction of the cubic perovskite-type reference structure. Our polarization analysis 
shows that partial contributions to polarization from the Fe and O atoms almost cancel each other 
and the net polarization present in BiFeOs mainly (> 98%) originates from Bi atoms. We found 
that the large scatter in experimentally reported polarization values in BiFeOs is associated with 
the large anisotropy in the spontaneous polarization. 

PACS numbers: 70., 78.20.Ls, 78.20.Ci, 74.25. Gz 



I. INTRODUCTION 



Multiferroic materials have coupled electric, magnetic 
and/or structural order parameters that result in simul- 
taneous ferromagnetic, ferroelectric and/or ferroelastic 
behavior. The magnetoelectric materials have two or- 
der parameters, viz. spontaneous polarization for ferro- 
electric or antiferroelectric cases and spontaneous mag- 
netization for ferromagnetic, ferrimagnetic or antiferro- 
magnetic cases. Coupling between the spontaneous po- 
larization and spontaneous magnetization lead to mag- 
netoelectric (ME) effects, in which the magnetization 
can be tuned by an applied electric field and vice versa. 
There has been considerable recent interest [11, [2, Sj in 
developing multiferroic materials in which interaction be- 
tween the magnetic and electronic degrees of freedom 
allows both charge and spin to be manipulated by ap- 
plied electric/magnetic field [1, These materials have 
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opportunities for potential applications in information 
storage, the emerging field of spintronics, and sensors. 
Other potential applications include multiple state mem- 
ory elements, in which data is stored both in the electric 
and the magnetic polarizations, or novel memory me- 
dia which might allow writing of a ferroelectric data bit, 
and reading of the magnetic field generated by associa- 
tion. Besides the potential applications, the fundamental 
physics of magnetoelectric materials is rich and fascinat- 
ing. Hence, the theoretical understanding of the intrin- 
sic physical properties of multiferroic oxides is clearly of 
great importance for both fundamental science and tech- 
nological applications. 

BiFeOs is ferroelectric with relatively high Curie tem- 
perature (Tc]ca. 1100 K) and exhibits antiferromag- 
netic behavior with high Neel temperature {T^;ca. 643 
The Fe magnetic moments are coupled fer- 
romagnetically (F) within the pseudocubic (111) planes 
and antiferromagnetically (AF) between adjacent planes. 
If the magnetic moments are oriented perpendicular to 
the [111] direction, the symmetry also permits a canting 
of the AF sublattices resulting in a macroscopic mag- 
netization; so-called weak ferromagnetism [13] ■ The 
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magnetic structure of BiFeOa is, in the first approxima- 
tion, AF with G-type magnetic ordering (G-AF). But the 
G-AF has been modified by subjecting it to a long range- 
modulation, manifesting itself in a cycloidal spiral of A= 
620 A length with [110] as the spiral propagation direc- 
tion and spin rotation within (110) which is unusual for 
perovskites [ll[. BiFeOa shows the linear magnetoelec- 
tric effect, with applied magnetic fields inducing weak 
F and large increase in polarization [l3|- The exhibi- 
tion of weak magnetism at room temperature is due to 
a residual moment from the canted spin structure Q • A 
significant magnetization {~ l/x^ per unit cell), as well 
as a strong ME coupling, have been reported recently in 
high quality thin films [13J. Coexistence of ferroelectric- 
ity and ferromagnetism has been observed [l^ in bulk 
BiFeOa-PrFeOa-PbTiOs solid solutions with composi- 
tions 0.2BiFeO3-0.2PrFeO3-0.6PbTiO3 and 0.4BiFeO3- 
0.2PrFeO3-0.6PbTiO3. Significant dielectric enhance- 
ment in ferrornagnetic 0.3BiFcO3-0.7SrBi2Nb2O9 has 
been reported [l^, for samples synthesized by sinter- 
ing mechanically activated oxide mixture. Since epitax- 
ial films of BiFe03 also show large electric polarization 
(~ 158/^ C/cm2) 16^ they are promising candidate mate- 
rials for ME-device applications. 

An important aspect that 
of the 
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BiFe03 is that it has a spontaneous polarization that 
is significantly smaller than the expected value for a 
ferroelectric with such a high Tc- It was until recently 
not clear whether this is a result of intrinsic material 
properties or of limitations imposed by leakage and 
imperfect material quality in bulk. Enhancement of 
spontaneous polarization and related properties have 
been reported [H, [H, [ll| for heteroepitaxially con- 
strained thin films of BiFeOs recently and more recent 
theoretical calculations [l^, llfll also predicated correctly 
the presence of large polarization in this material. The 
sizable polarizations in BiFe03 are consistent with 
the observed large atomic distortions [2l|, but 
apparently inconsistent with earlier studies of bulk 
BiFe03 SI, HI, origin of these distinction being 
currently under debate. This is one of the motivations 
for the present study. 

Recent theoretical report on structural behaviors of 
BiFe03 using density functional theory (DFT) calcula- 
tions show considerable deviation from the experimental 
findings [H, [l^ . An ideal cubic perovskite-type struc- 
ture can alternatively be described in a rhombohedral 
cell with rhombohedral angle {au)= 60°. The low tem- 
perature neutron measurements ^25i] show = 59.344° 
and 59.35° from single crystal X-ray diffraction measure- 
ments [l^l- This can be rationalized as a elongation of 
an ideal cubic perovskite-type structure along the [111] 
direction. In contrast, the hitherto available theoretical 
results available [l^ [l9| for BiFe03 conclude with an = 
60.36° which requires compression along the [111] direc- 
tion of the cubic perovskite like lattice, just opposite to 
the experimental findings. In order to find out whether 



DFT has failed to reproduce the structural parameters 
correctly for BiFe03, accurate calculations are needed. 
This is another motivation for pursuing this work. 

Moreover, observation of enhanced polarization in het- 
eroepitaxially constrained thin films of BiFe03, and the 
structural analysis [l^ suggest an approximately tetrag- 
onal structural arrangement with a superimposed mono- 
clinic distortion of about 0.5°. So, it is indeed interesting 
to look for possible metastable phases with slightly dis- 
torted structural arrangements to identify stable struc- 
ture of the experimentally-obtained heteroepitaxially- 
constrained thin films. Also, structural stability studies 
on BiFe03 would be useful to identify the paraelectric 
phase. Total energy studies on ferroelectric BiFe03 in 
different magnetic configurations would give more insight 
into its magnetic properties. Hence, we have made full 
structural optimization for BiFe03 in different atomic 
arrangements also including different magnetic configu- 
rations in the calculations. 

The "lone pair" in Bi or Pb-based oxides is believed to 
form due to the hybridization of 6s and 6p atomic orbitals 
with 6s^ electrons filling one of the resulting orbitals in 
Bi or Pb in their oxides. The lone pair is then considered 
to be chemically inactive, not taking part in the forma- 
tion of bonds but sterically active. The hybridization 
causes the lone pair to lose its spherical symmetry and 
is projected out on one side of the cation, resulting in an 
asymmetry of the metal coordination and distorted crys- 
tal structures. In the present study we have analyzed the 
occurrence of lone pair and its role on structural stability 
with the help of charge density, charge transfer, electron 
localization function, crystal orbital Hamiltonian popu- 
lation, and partial density of states analyses. 

The long-range Coulomb interaction plays a crucial 
role in ferroelectric materials. It is therefore critical to 
be able to calculate the spontaneous (i.e., zero field) elec- 
tric polarization P as a function of structural degrees of 
freedom. It has been shown [2^, [23] that P can directly 
be calculated as a Berry phase of the Bloch states. As 
for the dependence of P upon the structure, the first- 
order variation of P with atomic displacements is given 
by the dynamic effective charge Z*, which can be com- 
puted either directly using linear-response methods [2^ 
or numerically by finite differences. For the present study 
Born effective charges were obtained from finite differ- 
ences of macroscopic polarization induced by small dis- 
placements of the atomic sublattices. Recent theoretical 
reports [l^, [2^ indicated that it is important to include 
Coulomb correlation effects via the LDA+U approach 
into the calculation to describe the electronic structure 
and ferroelectric behavior of BiFe03. This motivated us 
to perform accurate DFT calculation to check the possi- 
bility of describing electronic structure, magnetism, and 
ferroelectric behavior of BiFe03 within the usual DFT 
framework. 

In this paper we report results of first principle DFT 
calculations on structural optimization for BiFe03 in dif- 
ferent atomic arrangements and also with different mag- 
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netic configurations. Moreover, we investigate structural 
phase stability, the electronic structure, chemical bond- 
ing, magnetism, and ferroelectric properties of BiFeOs. 
The remainder of this paper is organized as follows. In 
section |lTl we describe the structural aspects and the 
computations involved in the evaluation of structural 
phase stability, magnetism, and ferroelectric behaviors. 
In Sec. mil we present our results on structural phase 
stability, electronic structure, chemical bonding, mag- 
netism, and ferroelectricity, and compare the findings 
with available experimental and theoretical studies. We 
summarize and conclude in Sec. lIII II 



II. COMPUTATIONAL DETAILS 

First principles DFT calculations were performed us- 
ing the Vienna abinitio simulation package (VASP) 3(1 
within the projector augmented wave (PAW) method '31l| 
as implemented by Kresse and Joubert js^]. The Kohn- 
Sham equations [33| were solved self-consistently using 
an iterative matrix diagonalization method. This is 
based on a band-by-band preconditioned conjugate gra- 
dient [sj] method with an improved Pulay mixing |35j | 
to efficiently obtain the ground-state electronic struc- 
ture. The forces on the atoms were calculated using the 
Hellmann-Feynman theorem and they are used to per- 
form a conjugate gradient relaxation. Structural opti- 
mizations were continued until the forces on the atoms 
had converged to less than 1 meV/A and the pressure 
on the cell had minimized within the constraint of con- 
stant volume. The calculations were performed within 
periodic boundary conditions allowing the expansion of 
the crystal wave functions in terms of a plane- wave basis 
set. 

The Generalized Gradient Approximation (GGA) [s^ 
includes the effects of local gradients in the charge density 
for each point which generally gives better equilibrium 
structural parameters than the local density approxima- 
tion (LDA). Ferroelectric properties are extremely sen- 
sitive to structural parameters (lattice parameters and 
atom positions) and hence we have used GGA for all 
our studies. In the basis we treated explicitly 15 va- 
lence electrons for Bi {5d^°6s^6p^), 14 for Fe (3/3d^4s^ 
and 6 for oxygen {2s'^2p'^). Brillouin zone integrations 
are performed with a Gaussian broadening [33] of 0.1 eV 
during all structural optimizations. These calculations 
are performed with a 4x4x4 Monkhorst-Pack k-point 
meshjss'l centered at F for the ferroelectric i?3c struc- 
ture and similar k-point density has been used for all the 
other structures considered for the present study. One 
specific problem of these materials is the presence of 
"computationally-unfriendly" atom configurations such 
as the first-row element oxygen and the transition metal 
iron. Both are kinds of atoms which require extra care: 
either large basis sets within a pseudopotential scheme, 
or an all-electron scheme. So we have used very large 
basis set with 875 eV for the plane-wave cutoff. For the 



k-space integrations in the Berry-phase calculations, a 
uniform 8x8x8 k-point mesh was found to be ade- 
quate. 

BiFeOs has a rhombohedrally distorted perovskite- 
type structure with space group R3c [111, For 
the study of structural phase stability we have consid- 
ered eight closely related potential structure "types" : 
BiFeOa (rhombohedral; R3 c] | 23 | ■ CaTiOa (perovskite 
prototype; cubic;Fm3w) |3g |. KNbOa (tetragonal; 
PAmm) /3-LaCo03 (rhombohedral; i?3c) [il[, a- 

LaCoOa (monoclinic;/2/a) [i^, PbFeo.5Nbo.5O3 (mon- 
oclinic; Cm) [4^, GdFeOa (orthorhombic Pnma), and 
LaVOa (monoclinic; P2i/c) [4j. The reasons for choos- 
ing these structure types for our calculations are given 
below. It is experimentally reported that BiFeOa sta- 
bilizes in the rhombohedral i?3c and hence it is natural 
to include this structure variant. This structure which is 
also adopted by the paraelectric /J-LaCoOa structure can 
be derived from the ideal cubic perovskite-type structure 
by rotation of the octahedra along the cubic [111] axis. 
Since the experimental studies show enhancement of the 
polarization in heteroepitaxially constrained thin films of 
BiFeOa , for which the structural analyses [11] suggest an 
approximately tetragonal crystal structure with a small 
monoclinic distortion, we considered structural arrange- 
ments within space groups PAmm and Cm. Recent ex- 
perimental studies show that the /J-LaCoOa-type rhom- 
bohedral structure undergoes a small monoclinic distor- 
tion and stabilizes in an a-LaCo03 type structure [4^ . 
The orthorhombic GdFeOa-type structure arises from the 
ideal cubic perovskite by tilting of the FeOe octahedra 
around the cubic [110] axis and with small monoclinic 
distortion one arrives at the LaVOa-type structure [i^ . 
As experimentally reported crystal structure of thin film 
has small monoclinic distortion, we have considered the 
monoclinic I2/a, Cm, and P2i/c structures in our inves- 
tigations. 

The structures were fully relaxed for all volumes con- 
sidered in the present calculations using force as well as 
stress minimization. Experimentally established struc- 
tural data were used as input for the calculations. In or- 
der to avoid ambiguities regarding the free-energy results 
we have always used the same energy cutoff and corre- 
sponding k-grid densities for convergence in all calcula- 
tions. A criterion of at least 0.01 meV/atom was placed 
on the self-consistent convergence of the total energy. 
Brillouin zone integration was performed with a Gaus- 
sian broadening of 0.1 eV during all relaxations. The 
present type of theoretical approach has recently been 
successfully applied [i^, [4^ to reproduce experimentally 
observed ambient- and high-pressure phases of metal hy- 
drides. 

The experimental studies show that the ground state of 
BiFeOa is AF with a Neel temperature of 643 K [13, 113 ■ 
A noticeable magnetization (~l/iB per cell), as well as a 
strong magneto-electric coupling have been reported re- 
cently in high-quality epitaxial thin films [isj . The Fe 
magnetic moments are coupled F within the pseudocu- 
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bic (111) planes and AF between adjacent planes. It 
has been reported as a long-wavelength spiral spin struc- 
ture [m and possibly a small out-of-plane canting due 
to weak ferromagnetism [4^. Owing to tilting of FeOe 
octahedra, Fe-O-Fe chains are not subtending to an an- 
gle of 180°. As the exchange interaction takes place 
via these Fe-O-Fe chains, canted spin arrangements arise 
which further lead to antiferromagnetism with weak fer- 
romagnetism. Since the noncoUinearity of the magnetic 
moments appears to be quite minimal, we made the com- 
putational simplification to approximate the magnetic 
structure to a coUinear model. 

BiFeOs has three possible AF arrangements depending 
on the interplane and intraplane couplings within struc- 
ture, (i) With interplane AF coupling and intraplane F 
coupling the A-AF structure arises, (ii) The opposite 
structure, with the interplane F coupling and intraplane 
AF coupling is called C-AF arrangement. In the C-AF 
magnetic structure all magnetic atoms have two F and 
four AF nearest neighbors (vice-versa in the case of A- 
AF). (iii) If both the inter- and intraplane couplings are 
AF the G-AF arrangement arises. In a G-AF arrange- 
ment for an ideal cubic perovskite-type structure, each 
magnetic atom is surrounded by six AF neighbors. Be- 
cause of the rhombohedral distortion we were unable to 
introduce perfect G-AF ordering (see Fig. [T|). However, 
the adopted G-AF magnetic arrangement has perfect 
AF ordering according to cubic perovskite-type structure 
along the ac plane and zig-zag chains with AF ordering 
in the be plane. For the structural phase stability calcula- 
tions we assumed AF ordering with q vector (0,0,1/2) for 
all structures except for the cubic variant where we have 
assumed G-AF ordering and for the R3c variant where 
we considered nonmagnetic, F, A-AF, C-AF and G-AF 
orderings. 

In order to understand the role of Bi lone pair on the 
ferroelectric properties of BiFeOa we have visualized it 
with the help of electron localization function (ELF). The 
ELF provides a measure of the local influence of the Pauli 
repulsion on the behavior of the electrons and permits 
mapping in real space of core, bonding, and nonbonding 
regions in a crystal. It can be used to distinguish the 
nature of different types of bonding in solids [43]. It is 
defined as 

ELF=[l + i^r]-\ (1) 

where 

i 

and 

D, = A(3vr2)5/3p5/3 

here p is the electron density and (f)i are the Kohn-Sham 
wave functions. 



III. RESULTS AND DISCUSSIONS 
A. Structural phase stability 

BiFeOa is reported to exhibit about eight struc- 
tural transitions and a weak F ordering in its ground 
state structure [13, HH- The structural instability in 
perovskite-like oxides is partly related to the lattice mis- 
match between the cations. Such a mismatch is conve- 
niently quantified in terms of the empirical Goldsmith 
tolerance factor t. It is defined as the ratio between the 
ideal cubic lattice parameters based on A— O vs B— O 
bonding alone, with t = 1 corresponding to high sta- 
bility of the cubic perovskite structure. Using the ionic 
radii given by Shannon and Prewitt 52], the calculated 
tolerance factor [t — {tbi + tq /^/2{rFe + f o)] for BiFeOa 
is 0.96. Relatively low value of t for BiFeOa suggests 
that Bi ion may be unstable in the high-symmetric lo- 
cation possessed by the cubic perovskite-type structure, 
viz. consistent with the total energy results discussed 
below(shown in Fig:[2). In order to minimize the lattice 
mismatch and gain stability, a lattice distortion arises 
by the cooperative rotation of FeOg octahedra along the 
Jill] direction. In the ideal cubic perovskite structure, 
the BiOa planes are equidistant from the neighboring Fe 
(111) planes. Due to the rotation of the octahedra, three 
Bi— O distances are reduced, and three increased, giv- 
ing rise to an octahedral distortion with rhombohedral 
symmetry. 

The structure of the ferroelectric BiFeOs phase has 
been resolved experimentally using x-ray and neutron 
diffraction [13, HJ, [H, S?!, [H, [H, , and found to pos- 
sess a highly distorted perovskite structure with rhombo- 
hedral symmetry and space group i?3c. The unit cell of 
BiFeOa arising from the cubic perovskite structure by 
two kinds of distortions is shown in Fig.[3l One dis- 
tortion is the polar displacements of all the anion and 
cation sublattices relative to each other, which lead to 
the spontaneous electric polarization, and the other is 
an antiferro-distortive rotation of the FeOg octahedra 
along the Jill] direction with alternating sense of rota- 
tion along the Jill] axis (owing to the lattice mismatch 
as mentioned above). In terms of symmetry groups, the 
polar displacements alone would reduce the symmetry 
from cubic Pm3m to rhombohedral R3m, whereas the 
rotation of the FeOg octahedra alone would lead to the 
paraelectric phase with the space group i?3c. The in- 
corporation of both kinds of distortions gives the actual 
ferroelectric phase of BiFeOs with the space group R3c. 

The calculated total energy vs volume for BiFeOs in 
fully relaxed potential structures in the AF states are 
given in Fig.O This figure shows that the experimentally 
observed ferroelectric i?3c phase represents the ground 
state structure for BiFeOa. The calculated equilibrium 
structural parameters for BiFeOs in the ground state and 
closely related structures are given in Table|Tl Our cal- 
culated equilibrium structural parameters are also com- 
pared with available experimental and theoretical values 
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FIG. 1: (Color online) Antiferromagnetic structures considered for ferroelectric BiFeOa in the (a) A-AF-like structure, (b) 
C-AF-like structure , and (c) G-AF-like structure. For clarity only the Fe atoms are shown in a 2 x 2 x 2 of the primitive cell 
used in the calculation. 
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FIG. 2: Calculated cell volume vs total energy for antiferro- 
magnetic BiFeOs in different possible structural arrangements 
(structure types being labelled on the illustration). The ar- 
row indicates the pressure-induced structural transition point 
from the rhombohedral (_R3c) to the orthorhombic (Pnma) 
structure. 



in this table. The equihbrium volume obtained from the 
calculation is overestimated compared to the experimen- 
tal volume by ~ 3% and the lattice parameters are overes- 
timated by only 1.1%. The overestimation of equilibrium 
volume by GGA for rhombohderal phase of ferroelectric 
materials has already been observed for BaTiOS and Pb- 
Ti03 (56j . small deviation between experimental and the- 
oretical values indicate the predictive capabilities of the 
present type of state-of-the-art density functional calcu- 
lations. It should be noted that the same method with 
lower energy cut-off and k-point densities yielded the 




FIG. 3: (Color online) Crystal structure of BiFeOa in the fer- 
roelectric R3c structures. Highly distorted FeOe octahedra 
are corner shared through oxygen. Also note that Bi atoms 
are not present exactly in the middle between the FeOs oc- 
tahedra owing to the off-center displacement of Bi due to the 
presence of 6s lone pair electrons. 



equilibrium volume which is much smaller (~ 6.9%) than 
the experimental values. An earlier report [13,] shows that 
the ferroelectric properties of BiFeOs are very sensitive to 
the small changes in the lattice parameters. So, accurate 
calculations for the estimation of structural parameters 
for BiFeOa are needed in order to predict the ferroelectric 
properties correctly. 

The hypothetical nonferroelectric cubic perovskite- 
type structure is found to exhibit the highest energy 
among the various BiFeOa phases tested in this study. 
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The energy difference between tlie prototypical cubic 
structure and tlie equilibrium ferroelectric R3c structure 
is around 1.1 eV/f.u. This indicates that the paraelectric 
phase of BiFcOa does not take the cubic perovskite-like 
structure. In fact several low symmetric atomic arrange- 
ments based on tetragonal, rhombohedral, and mono- 
clinic structures are found at lower total energy than the 
cubic perovskite-type structure (see Fig.[2|). It is wor- 
thy to note that several low symmetry structures are 
found to be stable at high pressures in the closely re- 
lated system PbTiOa [s^l ■ The stabilization of low sym- 
metric structures is associated with the presence of lone 
pair electrons on Bi and also the hybridization between 
Bi and O. At expanded lattice we have identified two 
ferroelectric structures which are energetically closer to 
the ground state structure. One is a tetragonal PAmm 
structure and another is a monoclinic Cm structure. Ac- 
cording to a recent structural refinement by Noheda ct 
al. nil PbZro.52Tio.48O3 (PZT) at 20 K is in the mono- 
clinic Cm structure. PZT also stabilizes in the tetragonal 
PAmm structure. This material exhibits extremely large 
electromechanical coupling in both tetragonal and mono- 
clinic phases. It is also interesting to note that enhance- 
ment of spontaneous polarization and related properties 
in heteroepitaxially constrained thin films of BiPeOs has 
been recently reported [H, [l^, [is'l . In this context obser- 
vation of metastable Cm and PAmm phases in BiFeOs 
is particularly interesting. 

The total energy curves in Fig.H] show that there is 
a large difference in total energy between considered 
phases at the expanded lattice. However, at high pres- 
sures the energy difference is small which is associated 
with the weakening of directional bonding at higher 
pressures. Interestingly we observed a pressure-induced 
structural transition from the ground state rhombohedral 
R3c structure to a monoclinic P2i/c structure around 
13 GPa. As the monoclinic distortion was negligibly 
small, our detailed analysis of the high pressure phase 
yielded an orthorhombic Pnma (GdFeOa-type) struc- 
ture. The predicted structural parameters for the high 
pressure phase at the phase transition point are given in 
Table H 

The calculated bulk modulus for the ground state 
phase and the possible meta stable phases are also given 
in Table HI The bulk moduli have been obtained using 
the so-called universal equation of state fit for the total 
energy as a function of volume. Experimental bulk mod- 
ulus values are not available for this compound. The 
calculated bulk modulus for the ground state and the 
high pressure phases are found to be almost double the 
value of that for the other metastable phases given in Ta- 
ble U and this is associated with the smaller equilibrium 
volume of the former phases. 

The ideal cubic perovskite structure can be consid- 
ered as a special case of the rhombohedral structure with 
a rhombohedral angle of 60°. The experimental stud- 
ies [H Hi] show that an for BiFeOa is 59.344°. In con- 
trast, the hitherto available theoretical results [l^ [l^ 
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FIG. 4: Total energy as a function of rhombohedral angle 
(qh) for BiFeOs in the ferroelectric i?3c phase. 



show that the an value is 60.36° for BiFeOa which can 
only be achieved by compression of the cubic perovskite- 
like lattice along the [111] direction. With this back- 
ground it is interesting to calculate the variation in to- 
tal energy as a function of the rhombohedral angle for 
BiPeOa for the equilibrium volume. The calculated total 
energy vs an given in Fig. [3] clearly show that the the- 
oretical equilibrium rhombohedral angle obtained from 
accurate total energy calculation is much smaller than 
60° which is consistent with experimental observation. 
Compared with the total er ierg y corresponding to the ear- 
lier theoretically reported [ij, [l^l value for ur, a, gain of 
12.9meV/f.u. is obtained (the total energy gain from 60° 
to 59.235° is 6.5 meV/f.u.). Hence, our theoretical study 
confirms the experimental observation that the rhombo- 
hedral distortion in the ferroelectric BiFeOs phase cor- 
responds to an elongation rather than a compression of 
cubic perovskite-type lattice along [111]. 

We have also measured high resolution synchrotron 
powder diffraction data collected at the Swiss- Norwegian 
beam lines at ESRF, Grenoble, for the 29 range from 3 
— 38° and A= 0.49969 A. The structural parameters ex- 
perimentally measured at 115 K are given in TableUand 
found to be in good agreement with our calculated equi- 
librium structural parameters for the ferroelectric phase. 
We also found no significant difference in the measured 
structural parameters at 115 K for a sample subjected to 
0.8 T magnetic field during measurements. More details 
about the sample preparation and other experimental de- 
tails will be published later. 

In order to identify the possible polarization path we 
mapped the total energy as a function of displacement of 
Bi ions with respect to the FeOe octahedra along [001], 
[110] and [111] directions in the paraelectric R3c phase in 
Fig. [5] This helps to visualize the directions of the "eas- 
iest" ferroelectric transformation path and to clarify the 
role of Bi ions on the ferroelectric distortion. The poten- 
tial energy surface associated with the displacement of Bi 
ion along all the three directions (see Fig.E^) are having 
a double well shape. From Fig.[5]it is clear that the lowest 
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TABLE I: Optimized structural parameters and bulk modulus (Bo) for BiFeOa in different structural arrangements. 



Modification 


Unit cell 




Positional parameters 


Bo (GPa) 


(structure type) 


dimension (A) 


or ° 








-R3c (Present - theory) 


a = 5.697; a = 


- 59.235 


Bi(2a): 


0, 0, 0; Fe(2a): .2232, .2232, .2232 


130.9 




Vo = 128.48 




0(6b): 


.5342, .9357, .3865 




R3c (Ref. 19 - theory) 


a — 5.46; a — 


60.36 


Bi(2a): 


0, 0, 0; Fe(2a): .231, .231, .231 




Vo = 115.98 




0(6b): 


.542, .943, .408 




R3c (Ref. 59 - theory) 


a — 5.50; a — 


59.99 


Bi(2a): 


0, 0, 0; Fe(2a): .228, .228, .228 




Vo = 117.86 




0(6b): 


.542, .942, .368 




-R3c (Ref. 22 - experiment) 


a — 5.63; ct — 


59.35 


Bi(2a): 


0, 0, 0; Fe(2a): .221, .221, .221 




Vo = 124.60 




0(6b): 


.538, .933, .395 




R3c (Present - theory) 


a — 5.513; a - 


= 61.432 


Bi(2b): 


1/2, 1/2, 1/2; Fo(2a): 0, 0, 


89.8 




"\ 7 1 on nn 
V o — LZZ.Z)J 




0(6b): 


.3358, .1641 1/4 




R3c (Ref. 19 - theory) 


a — 5.35; a — 


61.93 


Bi(2b): 


1/2, 1/2, 1/2; Fc(2a): 0, 0, 




Vo = 113.12 




0(6b): 


.417, ..083 1/4 




Pnma (Present-theory) 


a = 5.5849; b= 


= 7.6597; c=5.3497 


Bi(4c): 


0.0536, 1/4, 0.98858; Fe(4b): 0, 0, 1/2 


138.8 




Vo = 113.12 




0(4c): 


.9750, 1/4, .4060; 0(8d): .2000, .9540 .1945 




PAmm (Present-theory) 


a = 3.7859; c= 


-- 4.8525; 


Bi(la): 


0, 0, .9491; Fo(lb): 1/2, 1/2, .51804 


72.3 




Vo = 139.10 




0(lb): 


1/2, 1/2, .1428; 0(2c): 0, 1/2, .6735 




Cm (Present-theory) 


a = 5.7900; b= 


-- 5.6899; c= 4.1739; 


Bi(2a): 


.9376, 0, .0685; Fo(2a): .5110, 0, ..4961 


79.0 




/3=91.99; Vo = 


137.42 


0(2a): 


.5626, 0, .9489; 0(4b): .7958, .7603, .4231 





energy off-center displacements are along the [111] direc- 
tions. This confirms that the polarizability of Bi plays 
a special role in the ferroelectric properties. The elec- 
tronic structure studies on Pb(Zro.5Tio.5)03 show (60j 
polarization rotation from the [100] to the [111] direc- 
tion owing to the smaller energy difference between the 
tetragonal and the rhombohedral structures. An experi- 
mental study [HH on PZT showed that the piezoelectric 
elongation in tetragonal PZT is along the direction asso- 
ciated with a monoclinic distortion, suggesting the mech- 
anism of polarization rotation via the monoclinic phase. 
Our total energy studies on BiFeOs show that the en- 
ergy difference between the Bi ion displacement along 
[111] and [110] directions are relatively small. Consistent 
with our observation, experimental studies indicate 
that heteroepitaxially strained thin films of BiFeOs have 
tetragonal-like crystal structure with the c axis normal to 
the substrate surface, with a small monoclinic distortion 
of about 0.5°. 

Figure[5]shows that the energy gain associated with the 
ferroelectric distortion along [111] direction is 0.427 eV 
per f.u. This energy gain is found to be much higher 
than 0.056 eV for PZT, the well known piezoelectric ma- 
terial [i^. It is tempting to try to connect these num- 
bers with the corresponding ferroelectric Curie temper- 
atures. We know that, in principle, the thermal energy 
required for a system to remain in the high-symmetric 
paraelectric i?3c phase will be of the order of energy bar- 
rier between these two systems. This suggests that the 
the Curie temperature of BiFeOa must be much higher 
than that of PZT. Consistent with this observation the 
ferroelectric Curie temperature of BiFeOa is (~ 1000 K) 
larger than that of PZT 570 K). It should be noted 
that for systems with a competition between different 
structural instabilities [11] , such an argument can grossly 
overestimate the transition temperature. In BiFeOa, for 
example, it is known that there is a competition between 
the ferroelectric rhombohedral phase (arising from the 
off-center displacements of ions due to Bi lone pairs) and 



an antiferro-distortive phase (associated with the lattice 
mismatch involving rotation of the oxygen octahedra). 

B. Chemical Bonding 

The ferroelectric properties of materials such as 
Pb(Zri_a;Tia;)03 are known to be related to the partial 
covalency of some bonds [6^ . [65j . In this connection de- 
tailed analysis of chemical bonding in BiFeOs is partic- 
ularly interesting. Analysis of bonding interaction from 
charge density distribution alone may mislead the identi- 
fication of the nature of chemical bonding correctly [6^ . 
Hence we have used different tools to explore the bond- 
ing behavior in BiFeOa. From the analyses of the charge 
density and densities of states for PbTiOs, and BaTiOa 
it is concluded [g^, [131 ^^at the ferroelectric instability 
is due to hybridization between the O 2p states and the 
Ti 3(i states. So, it is interesting to analyze the chemical 
bonding between Fe and O in BiFeOa in order to have 
a better understanding on the origin of ferroelectricity. 
The electron densities calculated by VASP were shown 
for a plane containing Bi, Fe and O atoms in Fig. [5^. 
From this figure it is clear that the bonding interaction 
between Fe and O is not purely ionic. The directional 
nature of charge density distribution between Fe and O 
indicates the presence of finite covalent bonding. 

Experimental results [11] and detailed theoretical in- 
vestigations [i^ on BaTiOs showed the significant hy- 
bridization interaction between Ba 5p and O 2p orbitals 
which contributes to the Born effective charge and fer- 
roelectric properties. In this context, the analysis of the 
nature of bonding interaction between Bi and O is par- 
ticularly interesting. The charge density distribution in 
FiglB^ shows presence of finite charges between Bi and 
O, which is generally an indication for the covalent inter- 
action. Charge transfer distribution (see FiglSJs) shows 
that electrons transfer from both Bi and Fe atoms to 
the oxygen sites consistent with the ionic picture. If one 
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Bi displacement (in A) Bi displacement (In A) 



FIG. 5: (a) Total energy as a function of displacement of Bi ion along [100], [110] and [111] directions in the R3c structure for 
the ferroelectric phase of BiFeOa. (b) The same for BiFeOa in the [111] direction with or without structural relaxation. 




FIG. 6: (Color online) Calculated (a) valence-electron charge density, (b)charge-transfer and (c) ELF plot for ferroelectric 
BiFeOs. 
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has pure ionic bonding interaction between Bi and O as 
well as between Fe and O, one would expect an isotropic 
charge transfer distribution. The anisotropic distribution 
of charge transfer confirms the presence of finite covalent 
bonding between Bi and O as well as between Fe and O. 

ELF is an informative tool to distinguish different 
bonding interactions in solids ^49] , here it is shown for fer- 
roelectric BiFeOa in Fig. [BJ;. The negligibly small value 
of ELF between atoms indicate the presence of domi- 
nant ionic bonding in this material. Also the ELF dis- 
tribution given in Fig. [6t shows maximum value at the 
O sites and minimum value at the Fe and Bi sites reit- 
erate charge transfer interaction from Bi/Fe to O sites. 
Moreover polarization of ELF at the O sites towards the 
other O sites and finite ELF between Bi and O indicate 
the hybridization interaction. The conclusion from the 
charge density, charge transfer and ELF analyses is that 
the bonding interaction between Bi and O as well as Fe 
and O are dominant ionic bonding with finite covalent 
character. This conclusion is consistent with the con- 
clusions arrived from the Born effective charge analysis 
discussed below. Our results clearly establish the mixed 
ionic-covalent character of the bonding in BiFeOs: it is 
obvious that the covalent character is not restricted to 
the Fe— O bond. It may be noted that while Bi s—0 p 
bonding is not very strong, it is at least as important as 
the Fe d—0 p interaction. 

If the bonding between the constituent atoms in 
BiFeOa is purely ionic it would remain centrosymmet- 
ric (and therefore not ferroelectric), because the short- 
range repulsion between adjacent closed-shell ions is min- 
imized for symmetric structures. The existence or ab- 
sence of ferroelectricity is determined by a balance be- 
tween these short-range repulsions, which favor the non- 
ferroelectric symmetric structure, and additional bond- 
ing considerations which act to stabilize the distortions 
necessary for the ferroelectric phase. Until now two dif- 
ferent chemical mechanisms have been identified for sta- 
bilizing the distorted structure in common ferroelectrics. 
In the first mechanism, the ferroelectric distortion was as- 
cribed to the hybridization interaction between the tran- 
sition metal ion and oxygen [6^ . The transition metal 
ions (such as Ti'*+, Zr^+, lSlh^+ and Ta^+) participat- 
ing in hybridization interaction with O are formally in a 
dP state, hence the lowest unoccupied energy levels are 
d-states which tend to hybridize with O 2p states [70| . 
However, this mechanism is not very significant for the 
case of multiferroics because the transition metal ions are 
not in a dP state. 

The second mechanism involves cations with "lone 
pair" electrons which have a formal ns^ valence electron 
configuration. In the same vein as the d*^ ions discussed 
above, these p° ions (T1+, Sn2+, Ph^+ , Sh^+ , and Bi^+) 
contain some p charge density which contribute to the 
displacive distortions. Indeed the tendency of ns^ ions 
to lose inversion symmetry is well established, with the 
conventional explanation invoking a mixing between the 
localized ns'^ states and a low-lying ns^np^ excited state, 




Energy (eV) 

FIG. 7: COHP for BiFeOa in the ferroelectric R3c structure 
describing Bi— O, Fe— O and Bi— Fe interactions. Ol and 02 
are the planar and apical oxygen atoms, respectively. 

a mixing that can onlyoccur if the ionic site does not have 
inversion symmetry [7ll | . If the lowering in energy asso- 
ciated with the hybridization interaction is larger than 
the inter- ionic repulsion opposing the ion shift, then a 
ferroelectric distortion results. This "stereochemical ac- 
tivity of the lone pair" is the driving force for off-center 
distortion in magneto-electric materials. 



C. COHP 

The crystal orbital Hamiltonian population (COHP) 
is calculated using the tight binding linear muffin-tin or- 
bital program (TBLMTO) The COHP is the density 
of states weighted by the corresponding Hamiltonian ma- 
trix element and is indicative of the strength and nature 
of a bonding (negative COHP) or antibonding (positive 
COHP) interaction 

The calculated COHP for the Bi-0, Bi-Fe, and Fe-0 
interaction for ferroelectric BiFeOa are shown in Fig.[7l 
The COHP curves reveal that all occupied states for the 
Fe— O interaction have bonding states, and the Fermi 
level is perfectly adjusted to fill up these bonding states. 
So the strongest bonding interaction is between Fe and 
O atoms. The FeOg octahedra are highly distorted and 
hence the Fe— 02 bonding energy (—1.52 eV) is stronger 
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than the Fe— 01 interaction (—1.02 eV) where 01 and 02 
are planar and apical oxygen atoms, respectively. Con- 
siderable amount of bonding states within the valence 
band (VB) in COHP shows significant bonding interac- 
tion between Bi and O atoms. It should be noted that a 
large bonding-state peak appearing around —10 eV (be- 
tween the Bi-0 interaction) is due to the hybridization 
interaction between the Bi 6s/6p and O 2s states. As the 
COHP is estimated from TBLMTO method, we found a 
shift in the DOS value by around 1 eV towards higher 
energy with respect to that from VASP. From the COHP 
analysis it is clear that the hybridization interaction be- 
tween Bi/Fe and O is significant, and it is an important 
driving force for the ferroelectric distortion. The COHP 
for Bi— Fc interaction shown in Fig. [7] indicates the pres- 
ence of both bonding and antibonding states within the 
VB. However, the bonding interaction between Bi and 
Fe in BiFeOa is weaker than that between all the other 
atoms. 



D. Band structure results 

In order to gain more insight into the role of ferro- 
electric distortion on electronic structure and bonding 
behavior in BiFeOs, the calculated band structure for A- 
AF phase of BiFeOs in the paraelectric i?3c structure and 
the ferroelectric R3c structure are given in Figs [5^ and 
b, respectively. As in R3c phase (Fig. [8^), the highest- 
energy valence states have noticeable O 2p character and 
the lowest conduction band states are Bi 6p-like, at vari- 
ance with the Ti 3d bands in PbTiOs. This together with 
the relatively narrow band gap of the paraelectric R3c 
phase, results in ferroelectric instabilities mainly char- 
acterized by a hybridization between Bi 6p and O 2p 
orbitals, not the one involving d bands of the transition 
metal cation that is typical of most other ferroelectric 
perovskites [H, It^I- It is interesting to note that the 
highest energy valence band state of i?3c phase has dom- 
inant Fe 3d states and these states are shifted to around 
0.35 eV to the lower energy in the paraelectric phase. In 
the R3c phase (see Fig.|S}D) the gap opens considerably 
and the top valence bands increase their Bi 6p character 
(See also Fig.[9l). The strong Bi 6p—0 2p hybridization 
is also reflected in the fact that the effective charge on Bi 
in the ferroelectric R3c phase has larger value than that 
in the paraelectric R3c phase (see Section llll Gl) . 

Two well dispersed bands around —10 eV is mainly 
originating from Bi-6s electrons. Compared with the 
paraelectric phase, these two bands are slightly broader 
in ferroelectric phase. The O 2s bands are completely 
filled and they are located around —18 eV below the 
Fermi level and are not shown in Figs. [5^ and b. The 
parabolic nature of bands around —6 eV are from Fe 4s 
electrons with significant contribution from both O 2p 
and Fe 3d electrons. A series of bands between —6 eV and 
the Fermi level are mainly originating from hybridized 
bands of both Fe 3c? and O 2p electrons. 



Owing to the off-center displacement of the ions in the 
ferroelectric phase of BiFeOs, the degeneracy in bands 
along the F — T direction in Fig. [5)3 are lifted compared 
to that of paraelectric phase. When the system is al- 
lowed to relax to its magnetic ground state a gap opens 
and an insulator is formed, as required for the ferroelec- 
tric state. The ferroelectric distortion slightly increases 
the band gap values, from 0.46 eV in the paraelectric to 
0.97 eV in the ferroelectric phases. The density func- 
tional calculations are known to underestimate the band 
gap and hence the reported values will be smaller than 
the real band gap value for this material. The Bi ions are 
almost in the 3-1- ionic state and hence the contribution 
from Bi 6p electrons are relatively small in the VB and 
they are mainly in the CB around 2 to 6 eV. 

In order to understand the sensitivity of electronic 
structure by ferroelectric displacement, the calculated to- 
tal density of states for the paraelectric and the ferroelec- 
tric phases of BiFeOs are superimposed in Fig. [S] From 
this figure it is clear that overall the energy of electrons 
in the ferroelectric phase are shifted around 0.5 eV from 
the paraelectric phase due to the off-center displacement 
of ions. This also highlights the different bonding inter- 
action between planar oxygen and Fe due to the buckling 
of the FeOg octahedra to accommodate the Bi 6s^ lone 
pair electrons occurring after the phase transition. In 
fact, the changes in the ferroelectric phase found at the 
top of the VB involve the d^z states of Fe and the Pz 
states of the apical oxygen, whose stronger hybridization 
gives more stability to the structure in the ferroelectric 
phase. 

To simplify the interpretation of the electronic struc- 
ture, we have also calculated the partial (atom and or- 
bital decomposed) densities of states for the ferroelec- 
tric BiFeOs (see Fig. [TU]). The top of VB is set to zero. 
Around —18 eV the lowest energy VB states form a filled, 
rather narrow band from O 2p electrons. Above these 
states (around —10 eV) lone pair Bi 6s states with small 
contribution from O 2s and 2p electrons are present. The 
Fe 4s states are located around —6 eV in the VB and they 
are completely mixed with the Fe 3d states in the entire 
VB. The energetically degenerate nature of Fe 3d and 
O 2p in the whole VB indicates the presence of covalent 
bonding between these two atoms. Bi 6s and 6p states 
are well separated and the latter states are mostly in the 
CB implying highly ionic nature of Bi ions. 



E. Magnetic properties 

Crystalline BiFeOa is a strong antiferromagnet (G-type 
spin structure, Tn = 654.9 K) [7a]- The cooperative mag- 
netism in BiFeOa is originating from the half-filled and 
localized (^igSg) Fe^+ ions. In an ideal cubic perovskite- 
like lattice, the Fe^^ ions in the high spin state prefer 
to form the G-AF ordering because, in this case the 
Pauli exclusion principle allows the transfer of electron 
to the neighboring ion in an antiparallel direction only. 
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FIG. 8; Electronic band dispersion for antiferromagnetic BiFeOs in (a) paraelectric R3c and (b) ferroelectric R3c structures. 
The Fermi level is set to zero. 
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FIG. 9: (Color online) Calculated total density of states for 
BiFeOa in the paraelectric R3c and ferroelectric R3c phases. 



The neutron diffraction studies carried out at room tem- 
perature revealed that BiFeOa, in fact, exhibits a com- 
pensated antiferromagnetic ordering (Tn = 670 K), with 
a cycloidal spin arrangement, incommensurate with its 
lattice [t^]. The magnetic structure of a separate unit 
cell corresponds to the G-type antiferromagnetic order- 
ing (each magnetic Fe'^"'" ion is surrounded by six Fe'^"'" 
ions with spins directed opposite to that of the central 
ion) 47]. Indirect evidence of the existence of the spa- 
tially spin-modulated structure (SSMS) in BiFeOa was 
obtained JUJ from the measurements of the magneto- 
electric effect. The existence of SSMS with periodicity 
of 620 A was recently confirmed by ^^Fe NMR measure- 
ments (78| . 

Dzyaloshinskii ^ developed a thermodynamic theory 
based on spin-lattice interaction and magnetic dipole in- 




-4 4 
Energy (eV) 

FIG. 10: (Color online) Partial density of states for ferroelec- 
tric BiFeOa in the G-AF state. 



teraction to explain the origin of weak ferromagnetism 
in antiferromagnetic compounds. The basic assumption 
is that antiferromagnetically ordered spins could rotate 
about a crystal axis toward one another, resulting in a net 
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TABLE II: Total energy (relative to the lowest energy state 
in meV/f.u.), Fe-magnetic moment (in HB/atom), and total 
magnetic moment (including those at oxygen and interstitial 
sites; /is/f.u.) for BiFeOa in A-, C-, G-AF, and F states in 
the ferroelectric R3c structure. 





A-AF 


C-AF 


G-AF 


F 


Energy 





211 


171 


390 


Fe-Magnetic moment 


3.724 


3.824 


3.789 


3.936 


Total magnetic moment 


3.885 


3.974 


4.076 


4.893 



spontaneous moment over the unit cell. The calculated 
magnetic moments and the total energy with respect to 
the ground state magnetic configuration for BiFeOs are 
given in Tablellll It is interesting to note that the en- 
ergy difi'erence between the ferromagnetic state and dif- 
ferent antiferromagnetic state is larger compared with 
that for LaMnOs (in LaMnOa the energy difference be- 
tween the ferromagnetic and the A-AF state is only 25 
meV/f.u.) [t^. This observation is consistent with the 
experimental measurements in the sense that amorphous 
phase of BiFeOa is found to be in a speromagnetic state 
(short range antiferro mag netic behavior) rather than in 
a ferromagnetic state [80| ■ 

Like LaFeOa, BiFeOs also has Fe ions in the 3-1- states. 
Similar to LaFeOa one could expect G-AF ordering in 
BiFeOa. But, owing to the presence of lone pair electrons 
at the Bi sites, one could expect anisotropy in the ex- 
change interaction also. Consistent with this expectation 
we found that the A-AF state is found to be 171 meV/f.u. 
lower in energy than the G-AF state in BiFeOa. This 
suggests that the Bi lone pair electrons play an impor- 
tant role in not only the ferroelectric behaviors but also 
in the magnetic properties of BiFeOa. The observation 
of A-AF ordering is particularly interesting. In LaMnOa 
each MnOe octahedron is elongated along one axis as a 
result of the Jahn— Teller distortion around the Mn 
ions. This Jahn- Teller distortion creates orbital order- 
ing that leads to four singly occupied d1 orbitals point- 
ing towards empty d^2_y2 orbitals on adjacent Mn ions 
(a ferromagnetic interaction) and two empty (1^2 _y2 or- 
bitals pointing towards other empty dx2-y2 orbitals (an 
antiferromagnetic interaction) resulting in the A-AF or- 
dering [roj . If no other mechanism was involved in the 
stabilization of magnetic ordering in BiFeOa one could 
expect G-AF ordering as that in LaFeOa. But, owing to 
the Bi lone pair electrons, we found A-AF phase to be the 
ground state for BiFeOa among the considered coUinear 
magnetic configurations. It should be noted that if the 
Coulomb correlation effects are important in this ma- 
terial, one could expect that the energetics of different 
magnetic configurations will be different from the values 
reported here. But, as most of the features of magneto- 
electric behavior of BiFeOS are well described by GGA 
itself, we have not considered Coulomb correlation effects 
in the present study. 

The calculated magnetic moment at the Fe site varies 
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FIG. 11: Calculated total density of states for BiFeOs in dif- 
ferent magnetic configurations. The Fermi level is set to zero. 



between 3.724—3.936 per Fe atom depending upon 
the magnetic configuration considered in our calculations 
(see Tableini. The calculated magnetic moment at the 
Fe site is found to be in good agreement with the value 
of 3.75(2)/iB measured from low temperature neutron 
diffraction measurements p5j . The calculated magnetic 
moments at the iron sites are not integer values, since 
the Fe electrons have a hybridization interaction with 
the neighboring O ions. Because of this hybridization 
interaction we found that around 0.18 magnetic mo- 
ment is induced at each O site in the ferromagnetic state 
which are polarized along the same direction as that in 
the Fe sites. As a result the net total moment at the site 
is considerably larger (see the total moment in Tablc ITT)) . 

The calculated total DOS for BiFeOa in the nonmag- 
netic and ferromagnetic states show metallic behavior 
with a large DOS at the Fermi level. The large num- 
ber of electrons near the Fermi level is not a favorable 
condition for stability. However, when one creates AF 
ordering, the exchange interaction produces an exchange 
potential that effectively shifts the energy of the Fe id 
band to lower energy giving insulating behavior as shown 
in Fig. 111! We found that even if we introduce spin 
polarization into the calculation, the electronic ground 
state depends upon the nature of magnetic ordering. The 
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Fermi level lies in a pseudogap feature in the C-AF state 
and clear insulating behavior is observed for the A-AF 
and G-AF ordered states with band gaps of 0.97 eV and 
0.53 eV, respectively. If the correlation effects are not sig- 
nificant in this material, the present result suggest that 
magnetic ordering plays an important role in stabiliza- 
tion of ferroelectricity. However, it should be noted that 
the magnetic ordering temperature is much smaller (643 
K) than the ferroelectric ordering temperature (1100 K). 
Hence, it appears that mechanisms other than magnetism 
may establish insulating behavior and ferroelectricity at 
high temperatures. 



F. Role of Bi lone pair electrons for ferroelectricity 

The ELF provides a measure for the local influence 
of the Pauli repulsion on the behavior of electrons and 
permits the mapping in real space of core, bonding, lone 
pair, and nonbonding regions in a crystal (8l| . We have 
used the TBLMTO program to calculate the ELF. In 
order to visualize the lone pairs around Bi, it was neces- 
sary to select an isosurface with an ELF value of 0.75 as 
shown in Fig. [T^l This figure reveals the lobe-like lone 
pairs possessed by the Bi atoms. It also indicates that 
the Bi 6s orbitals are in direct interaction with oxygen 2p 
and are thus not chemically inert. An interesting feature 
of ELF for BiFeOs in Fig. [3J; is the presence of a nonuni- 
form distribution of ELF at the Bi sites. The reason for 
this is that due to the displacement of the Bi atoms, the 
Bi 6pz orbital can interact more efhciently with the O 
2px and O 2py orbitals. In the R3c structure, these or- 
bitals interact only with the Bi 6s orbitals because the 
symmetry does not allow for an interaction with the Bi 
6pz orbitals. In the i?3c structure, half of the Bi atoms 
have moved in the +z direction and another half of them 
in the —z direction, breaking the symmetry and allowing 
hybridization between the indicated linear combination 
of O orbitals and the Bi p^ orbital. As a result, some of 
the O Px/0 Py derived states (in fact those that had the 
strongest interaction with the Bi 6s state) are lowered 
in energy and move away from the Fermi surface, effec- 
tively counterbalancing the destabilization caused by the 
Bi lone pair. The combination of large values of ELF and 
the covalency leads to a stronger lone pair, which is more 
stereochemically active. The manner in which the lone 
pairs dispose themselves suggests a greater tendency to 
the ferroelectric distortion. 

The lone pairs are expected to have primarily 6s char- 
acter. If they are to form lobes, there must be some 
hybridization interaction with the neighbors. The ques- 
tion is whether the p character arises from mixing s and 
p states on the Bi atom or is it an effect of covalency, aris- 
ing from Bi s and O s/p mixing. The partial DOS analy- 
sis (see Fig. 1101) shows that Bi 6s electrons are around 
— 10 eV below the Fermi level and are well localized. 
Also the Bi-s and p states are well separated indicating 
that the mixing between s and p on the Bi site is min- 




FIG. 12: (Color online) Isosurface (at a value of 0.75) of the 
valence electron localization function of BiFeOa in the ferro- 
electric i?3c structure obtained from TBLMTO calculation. 
The atom labels in the shown plane are same as that in Fig. [6^ 



imal. Watson and Parker 82] have recently suggested 
that it is due to p character from the anions that lone 
pairs are allowed to form lobe-shaped structures. The 
anion is clearly implicated since in a homologous series 
of compounds, changing the anion often greatly affects 
the stereochemistry of the lone pair [83| . If the lone pair 
6s electrons overlap with the O s/p states it can attain 
the typical lobe shape, conversely, in a situation where 
there is no such overlap, or when the overlap is weak, the 
lone pair could retain its pure 6s character and remain 
spherical. 

In BiFeOa the charge density and ELF analysis (see 
Fig. [6]) show noticeable covalent bonding between Bi and 
O. In order to have covalent bonding between Bi s and O 
p states they should be energetically degenerate. On the 
other hand the O p states are mainly present at the top 
of the VB. Although both Bi s and O s states are well 
localized and energetically well separated, they are spa- 
tially well dispersed. Owing to the spherically symmetric 
nature of the s orbitals, the probability of the hybridiza- 
tion interaction between s orbitals is much larger than 
that between any other orbitals. Hence the lobe shape 
of the lone pair is at least partly associated with the 
hybridization interaction between Bi s and O s states. 
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Even though O p states are mainly present at the top of 
the valence band, finite O p electrons are present around 
— lOeV where the Bi 6s electrons are also present. So, 
both O s and p electrons participate in the hybridization 
interaction with the Bi 6s electrons. Consistent with this 
view point our COHP analysis (see Fig. [7]) indicates that 
Bi— O bonding states are present around — lOeV. 

It was argued by Orgel [sJl, and has been assumed 
in much of the literature, that the interaction between 
occupied cation s and unoccupied cation p bands is the 
primary cause for the lobe formation. The partial DOS 
analysis clearly shows (see Fig. fTU]) that Bi s and p states 
are well separated. Hence our observation does not sup- 
port this classical view of hybridization of the Bi 6s and 
6p orbitals. Rather, it supports the observation of hy- 
bridization interaction between the anion s and cati on p 
states as seen in PbO and related compounds [H, [ssf . 
The formation of the distorted electron densities and 
crystal structures will thus depend on the energies of the 
valence orbitals of oxygen ions involved in bonding to the 
Bi. 

Although the s states of both Bi and O are well below 
the VB their importance for ferroelectric distortion can 
be understood as follows. The s states are spherically 
symmetric and diffuse compared to the p ov d states. 
Hence, the overlap interaction between s states will be 
stronger than that between other states. Owing to the 
highly ionic nature of Bi, the Bi Qp states are mostly in 
the unoccupied state (see Fig.[Tni) but with some notice- 
able Bi Qp states also at the VB. It has been shown that 
the stereochemical activity of the cation lone pair is de- 
termined by the interaction of cation s and p states with 
anion p states [SSj . Hence the activity of the lone pair will 
be stronger if the covalent interaction between Bi and O 
is larger, where it has the stereochemical effect of stabiliz- 
ing the distorted low symmetric structures. Stability of a 
distorted structure depends on whether the energy gained 
from hybridization interaction (bonding) between occu- 
pied and unoccupied states compensates the energy re- 
quired to promote an electron to unoccupied state. Since 
the same oxygen p states are involved in both interactions 
and hence in instabilities, a competition results. 

BiFeOs belongs to the ferroelectrics of displacive type. 
Now we try to understand the role of lone pair on the 
ferroelectric behavior. The lone pair occupies as much 
volume as an oxide or fluoride anion in crystals . The 
6s^ electrons of Bi^+ hybridize with 2s and 2p of oxygen 
to form a space- filling localized lobe, which in turn pushes 
away its neighboring atoms causing a structural distor- 
tion. The presence of lone pair makes the hybridization 
interaction between Bi and O different in different di- 
rections. Moreover, the effective ionic radius of Bi ions 
decreases on one side and increases on the other side. 
The anisotropic nature of electrons at the Bi sites (origi- 
nating from the presence of lone pair) stabilizes off-center 
displaced structure. This results in the movement of the 
Bi ion in the [111] direction which causes a cooperative 
displacement of the Fe'^"'' ions due to the repulsive inter- 



action and finally gives rise to a ferroelectric polarization. 

G. Born effective charge analysis 

The Born effective charge (BEG) Z* is a fundamental 
quantity for the study of lattice dynamics, because it gov- 
erns the amplitude of the long-range Coulomb interaction 
between nuclei, and the splitting between longitudinal 
optic (LO) and transverse optic (TO) phonon modes. In 
binary crystals ^87] , infra-red measurements of the split- 
ting between LO and TO modes allow an accurate esti- 
mation of ||Z*|p/eoo and offer therefore an unambiguous 
way to extract the amplitude of Z* from the experiment. 
But, in complex compounds like perovskite oxides, LO 
and TO mode eigenvectors are not necessarily equivalent 
and the determination of Z* from the experimental data 
is consequently not straightforward (88| . Hence, the de- 
velopment of theoretical methods giving direct access to 
Z* acquires specific interest. 

The dynamical charge measures the macroscopic cur- 
rent flowing across the sample while the ions are adiabati- 
cally displaced. Such currents are responsible for building 
up the spontaneous polarization, when ions are contin- 
uously displaced from the centrosymmetric structure to 
the ferroelectric structure 0, • These displacements 
are directly observable when measuring polarization via 
hysteresis cycles. The modern theory of macroscopic 
polarization [s^ allows the calculation of such currents 
as Berry phases. The dynamical Born effective-charge 
tensor Z* of a given ion measures (by definition) total 
macroscopic current adiabatically transported by the ion 
moving with unit velocity, while all the other ions and 
the macroscopic field are kept fixed (9(| . 

According to the modern theory of polarization [H, [H, 
[9l| , the change in macroscopic polarization between two 
different insulating states can be regarded as a measure of 
the phase difference between their initial and final many- 
body wavefunctions. The total polarization P for a given 
crystalline geometry can be calculated as the sum of ionic 
and electronic contributions. 

The ionic contribution is obtained by summing the 
product of the position of each ion in the unit cell with 
the nominal charge of its rigid ion core. The electronic 
contribution to P is determined by evaluating the phase 
of the product of overlaps between cell-periodic Bloch 
functions along a densely sampled string of neighboring 
points in k space. 

It has been shown [o^, [11] that the effective-charge 
tensors are almost constant along the path connecting 
the two structures, which implies that the spontaneous 
polarization is simply expressed in terms of the internal- 
strain displacement Ur as 



P, = (2) 
where VI is the cell volume, and Z* are the diagonal 
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components of the effective-charge tensors for atom r. 
BEC can be defined as a hnear term in the polarization 
change due to a unit displacement of the rth ion, keep- 
ing all the other ions fixed under the condition of zero 
macroscopic field i.e., 



^Tal3 — 



s=0 



(3) 



for the change in P on moving ion t; note that the effec- 
tive charge is a tensor. Thus the effective charge is cal- 
culated by finding the induced polarization of the nuclei 
and the electrons when the nuclei are displaced equally in 
each unit cell, so that the translational periodicity is the 
same as the undistorted crystal. The effective charge thus 
found is the transverse Born effective charge [o^, [o^ since 
the macroscopic electric field is required to vanish. Ac- 
cording to the modern theory of polarization [26,, ^89J , the 
total difference in polarization between a reference state 
and a state where the atoms belonging to the sublattice r 
have been displaced by a small but finite distance Aura 
is defined as, 



AP = AP, 



AP,:, 



(4) 



where APei is the electronic contribution obtained 
from the Berry-phase polarization approach and AP^on 
is the ionic contribution. APio„ is defined by 



AP,;, 



|e|zrAM 



(5) 



where Zr is the valence atomic number of the t*^ atom. 
For practical purposes, the electronic contribution to Z* 
can be obtained as 



lim 

AUr 



AP, 



+0 AUto 



(6) 



In periodic systems, the change in electronic polarization 
in zero field can be computed from the King-Smith and 
Vanderbilt formula 1261: 



'/3 ^ "7T^'^« '^'^ / {unk\ld/dkp\unk)dk (7) 
y^'^r Jbz 

where s is the occupation number of the valence band 
states (s = 1 for spin-polarized case) and Unk are the 
periodic parts of the Bloch functions. This definition is 
valid only under the constrain that the wave functions 
fulfil the periodic gauge condition. This means that the 
periodic parts of the Bloch functions must satisfy 



Unk{r) 



iG.r 
e Wnfe+G(r)- 



(8) 



Once AP is known, the Born effective charge tensor 
for the /? component can be obtained from the formula [3] 



where a denotes the direction of polarization. For a con- 
tinuous change, such as the induced polarization when 
atoms are displaced, uncertainties by possible integral 
multiples of 27r can be avoided by defining small enough 
changes where one always knows that the change in po- 
larization is given by the smallest value, i.e., with the in- 
teger equal zero. Hence, in the Berry-phase calculations, 
we have chosen the displacement of O.OSA, and used a set 
of strings of 8-k points (parallel to some chosen reciprocal 
lattice vectors) to calculate the electronic polarization. 
Our results shown in Table IIIII satisfy the acoustic sum 
rule, '^f, Zl ■■ = , indicating that the calculations are 
well converged with respect to computational conditions. 

Our calculated zz component of the BEC can directly 
be compared with the BEC values given in Table IV in 
Ref.[l^. The formal valence of Bi, Fe and O in BiFeOa are 
3-I-, 3-|- and 2— respectively. In simple, high-symm etry 
oxides, the oxygen Z* is isotropic and is close to —2 [9(| . 
However, simple models [97[ hint at nontrivial values of 
Z*. When we include AF ordering into the calculations, 
because of the loss of symmetry, the space group changed 
from R3c to i?3. Owing to the site symmetry, the diag- 
onal components of Z* are anisotropic for all the ions 
in Table imi If ions have closed-shell-like character, each 
ion will carry an effective charge close to their nominal 
ionic value according to a rigid-ion picture. On the con- 
trary, a large amount of nonrigid delocalized charge flows 
across this compound during the displacement of ions ow- 
ing to the covalence effect [63, US] • Consequently, we have 
obtained effective-charges much larger than the nominal 
ionic value. The charges on Bi, Fe and O are around 75% 
larger than they would have been in a pure ionic picture: 
this reveals the presence of a large dynamic contribution 
superimposed to the static charge. Similar giant values 
of Z* have been reported using quite different technical 
ingredients and/or for other perovskite oxides [69l. [ost. 

The BEC is indeed a macroscopic concept [13, UMl , in- 
volving the polarization of valence electrons as a whole, 
while the charge "belonging" to a given ion is an ill- 
defined concept. The high BEC values indicate that 
relative displacements of neighboring ions against each 
other trigger high polarization. Roughly speaking, a 
large amount of nonrigid, delocalized charge is respon- 
sible for higher value of BEC than the nominal charges. 
The calculated BEC for the ferroelectric and the para- 
electric phases are given in Table lllll It is clear that 
both Bi and Fe donate electrons and O accepts electrons, 
consistent with the traditional ionic picture. It can be re- 
called that for a pure ionic system one could expect more 
isotropic value for Born effective charges. Considerable 
anisotropy in the diagonal components of BEC given in 
Table HIIl and the noticeable off-diagonal components in 
the oxygen sites confirm the presence of covalent bonding 
between O 2p and Fe 3d orbitals. It is generally expected 
that there is considerable covalent bonding between Fe 
and O within the FeOg octahedra and this could explain 
the presence of anomalous contributions (defined as the 
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TABLE III; Calculated Born effective charge tensor for antiferromagnetic BiFeOa in the paraelectric R3c and ferroelectric R3c 



structures 



z* 


Position 


XX 


yy 


zz 


xy 


xz 


yx 


yz 


zx 


zy 


R3c 






















Bi 


(la) 


5.048 


4.975 


4.243 


0.104 


— 0.058 


—0.081 


— 0.058 


— 0.001 


0.000 


Bi 


(la) 


5.048 


5.071 


4.240 


— 0.105 


—0.058 


0.0330 


— 0.058 


0.000 


0.000 


Fe 


da) 


4.537 


4.607 


3.205 


— 0.185 


—0.069 


0.2170 


— 0.072 


0.000 


0.000 


Fe 


(la) 


4.537 


4.568 


3.207 


0.186 


-0.068 


-0.115 


-0.067 


0.000 


0.000 


O 


(3b) 


-3.799 


-2.601 


-2.583 


0.149 


0.265 


-0.038 


0.629 


0.211 


0.608 


O 


(3b) 


-2.978 


-3.469 


-2.584 


0.418 


0.434 


0.636 


0.554 


0.421 


0.486 


i?,3c 






















Bi 


(la) 


5.438 


5.460 


4.792 


-0.022 


-0.004 


-0.001 


0.004 


0.000 


0.000 


Bi 


(la) 


5.434 


5.414 


4.791 


0.022 


0.005 


0.001 


-0.005 


-0.001 


0.000 


Fc 


(la) 


5.224 


5.223 


4.859 


-0.777 


-0.001 


0.777 


0.000 


0.000 


0.000 


Fo 


(la) 


5.224 


5.222 


4.858 


0.777 


0.000 


-0.777 


0.000 


0.000 


0.000 


O 


(3b) 


-3.590 


-3.496 


-3.227 


-0.048 


0.409 


0.000 


0.000 


0.530 


-0.023 


O 


(3b) 


-3.551 


-3.533 


-3.226 


-0.062 


-0.212 


-0.032 


-0.347 


-0.285 


-0.446 



additional charge with respect to the well known ionic 
value) to the effective charge at the Fe and O sites. Sur- 
prisingly large anomalous charge at the Bi site (+1.8 to 
+2.4) gives tangible proof for hybridization interaction 
between Bi and the neighboring O atoms. When going 
from the paraelectric phase to the ferroelectric phase, 
we have observed considerable increase in the effective 
charges in all three constituents. This is associated with 
the anisotropy in the chemical environment. While FeOe 
octahedra are equally distanced from Bi ions in the para- 
electric phase, a short and a long Bi— FeOg distance arise 
in the ferroelectric phase owing to the off-center displace- 
ment of the Bi ions. 



H. Spontaneous polarization 

One of the basic quantities of ferroelectrics is spon- 
taneous macroscopic polarization P, which results upon 
application of an electric field, and persists at zero field 
in two (or more) enantiomorphous metastable states of 
the crystal. Experiments measure P via hysteresis cy- 
cles from the difference AP between these metastable 
states [9^. Even though BiFeOa has large atomic dis- 
placements relative to the centrosymmetric cubic per- 
ovskite structure, and the high ferroelectric Curie tem- 
perature (HOOK), early measurements yielded rather 
small polarizations [13, [H, H^l- The small values are 
in sharp contrast to t he ex perimental findings on epitax- 
ial thin films [l^ [l^ Il00l | which were found to possess 
large polarizations. There are several plausible explana- 
tions given for the large difference in values of polariza- 
tion reported experimentally. One possibility is that the 
original reports of small polarization might have been 
limited by poor sample quality [l7'|. A second possibil- 
ity is that the small values could be correct for the R'3c 
structure, with the large values being correct for different 
structural modifications stabilized in the thin films p^ . 
Recently the large and small values are explained within 
the modern theory of polarization, which recognized that 
polarization is in fact a lattice of values, rather than a 



single vector [19[. Hence it is interesting to analyze the 
spontaneous polarization present in BiFeOa in detail with 
the help of accurate density functional calculations. 

We calculated the polarization difference AP between 
the polar (i?3c; A=l) and the ideal (i?3c; A=0) rhombo- 
hedral structures assuming a continuous adiabatic trans- 
formation which consists of scaling the internal strain 
with the parameter A (0 < A < u). The polarization at 
the paraelectric R3c phase [P(0)] would be zero (since 
there is no dipole for this case) and would be a finite 
value P(u) for the ferroelectric i?3c phase. The dipole 
moment evolves continuously from P(O)=0 to P(u). i.e. 




Z*{u)du=[P{u)-P{0)]=uZ{u) (9) 



where 

1 r 

Z{u) ^ - / Z*{u)du (10) 
u Jo 

the mean value of Z*(u) from to u is equal to Z{u). 
We have calculated the polarization along x, y and z di- 
rections using three different approximations in Table HVl 
i.e. By assuming (i.) the Z{u) values same as Z*{u) 
values of ferroelectric phase, (ii.) the Z{u) values are 
same as Z*{u) values of paraelectric phase, and (in.) the 
Z{u) values is an average of Z*{u) value of the ferro-, 
and para-electric phases. The Z*^ eigenvalues of Bi and 
Fe correspond to an eigenvector aligned along the ferro- 
electric axis. Hence, the calculated polarization values 
corresponding to P^ obtained from (iii) can be directly 
compared with the experimentally-measured values. 

The i?3c symmetry permits the development of a spon- 
taneous polarization along [111], and Bi, Fe, and O are 
displaced relative to one another along this threefold 
axis. In the case of O, the eigenvector associated with 
the highest eigenvalue is positioned approximately along 
the direction of Fe— O bond. We identify that the high- 
est contribution originates from the equatorial plane of 
the FeOe octahedron. The largest relative displacements 
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TABLE IV: Cartesian components of the polarization A Pa 
for BiFeOs calculated using Eq[2] along a path from the para- 
electric R3c phase to the ferroelectric R3c phase. The polar- 
ization values (in /iC/cm'^) calculated with the Born effective 
charges for RSc, R3c structures and the average value for 
these two structures (noted by "average"). 



Z* values used in Eq.[2] 


P. 


P. 


Pz 


R3c 


-0.02 


-0.01 


93.85 


RSc 


1.33 


2.85 


83.53 


Average 


0.65 


1.41 


88.69 


Point charge 








60.19 



are those of Bi relative to O owin g to the presence of 
stereochemically active Bi lone pair |10l| . The off-center 
displacements are noticeably larger compared to those 
in nonlone-pair-activeperovskite ferroelectrics such as 
BaTiOa or LiNbOa [ij, Il02l |. but are consistent with 
those observed in other Bi-based perovskites [621] . 

Our partial polarization analysis shows that the polar- 
ization contribution coming from Bi, Fe, and O atoms 
are 86.94, 18.91, and — 17.15/iC/cm^ respectively along 
[111] direction. This indicates that the polarization con- 
tribution arising from displacement of Fe and O atoms 
are almost cancelling out and more than 98 % of the net 
polarization present in BiFeOs is contributed by the Bi 
ions. The oxygen ions are displaced in the same direction 
as Bi and Fe. However, the negative value of Born effec- 
tive charges at the oxygen sites is responsible for the op- 
posite direction of polarization at oxygen sites compared 
with that in Bi and Fe sites. We have also calculated the 
polarization using the point charge model and it is found 
to be around 32 % smaller than the actual value. The 
reduced value of polarization obtained from point charge 
model is related to the covalence effect. 

From Table HVl it is clear that the polarization val- 
ues are very much directional dependent (see Table ITV)). 
Our calculated value of polarization corresponding to the 
z direction can be compared with the theoretical val- 
ues of 84.2AiC/cm2 from LSDA [19] and 98.9 ^lC /cm^ 
from LDA+U method [1^. The large anisotropy in the 
calculated polarization values indicate that if the sam- 
ples are aligned along [111], then they will have larger 
polarization than those aligned along other orientations 
in the substrate. This could explain why some reports 
show large value (35.7 — 158/xC/cm^) of polarization 
in this material and other show very small value (2.2 — 
8.9/uC/cm^). We found that the natural value of the po- 
larization is 88.7/iC/cm^ along [111] direction, consistent 
with the recent measurements [13, Ug] . 



I. Summary 

We have reported results from electronic structure cal- 
culations for the multiferroic BiFeOa using generalized 



gradient corrected accurate density functional calcula- 
tion. We have found that A-AF-like phase is stable in 
energy compared to the other magnetic configurations in 
BiFeOa. This phase sustains ferroelectric rhombohedral 
distortions with large shifts of Bi atoms, in agreement 
with the fact that a large portion of the ferroelectric po- 
larization in Bi-based perovskite-like phases is provided 
by Bi displacements. The stereochemically active lone 
pair in BiFeOs causes a symmetry lowering structural 
distortion which drives the ferroelectric phase transition 
and, as a secondary effect, allows A-AF ordering or the 
canted spin structure over the G-AF ordering. The agree- 
ment between the calculated spontaneous polarization 
(P) and recent experimental data demonstrates that the 
crystal wave functions in the frozen ferroelectric structure 
convey the essential physics of the spontaneous polariza- 
tion in real materials. The important observations from 
our calculations can be summarized as follows. 

1. Unlike earlier calculations, our calculated structural 
parameters for ferroelectric BiFeOa are found to be in 
good agreement with experimental values and especially 
the experimentally-observed rhombohedral angle for this 
material has been reproduced quantitatively. 

2. Structural stability studies show a pressure induced 
structural transition from polar i?3c structure to the non- 
polar orthorhombic Pnma structure around 130kbar. 
We have also identified two ferroelectric phases such as 
tetragonal PAmm and monoclinic Cm structures at the 
expanded lattice. 

3. We have observed insulating behavior in both the 
paraelectric i?3c phase and the ferroelectric i?3c phase 
without including strong Coulomb correlation effect into 
the calculations. 

4. Total energy calculations show large energy difference 
between different magnetic configurations in BiFeOs and 
we found insulating behavior in A-AF and G-AF config- 
uration and metallic behavior in the nonmagnetic, fer- 
romagnetic and C-AF configurations. We believe that 
GGA itself sufficient to describe the properties of BiFeOs. 
However, if correlation effect is significant in this material 
one could expect that the predicted energetics of differ- 
ent magnetic configurations and electronic behavior will 
change. 

5. Detailed analysis of the chemical bonding shows that 
the electrons transfer from both Bi and Fe atoms to the O 
sites and also finite hybridization interaction takes place 
between the Bi-0 as well as Fe-0. This hybridization 
interaction plays a dominant role in the large ferroelectric 
polarization. 

6. The role of Bi lone pair on the chemical bonding and 
ferroelectric distortion has been analyzed. We found that 
the lone pair lobe formation is driven primarily by cova- 
lency between Bi s and O s/p states with the mediation of 
Bi p states rather than the classical picture of hybridiza- 
tion between Bi 6s and Bi 6p orbitals. 

7. The complete BEG tensor is reported for all the con- 
stituents in both the ferroelectric and paraelectric phases 
using accurate density functional calculations. The cal- 
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culated BEC at the ferroelectric and the paraelectric 
phase indicate that the charges are strongly affected by 
ferroelectric displacement of ions. 

8. Consistent with experimental observation, our calcu- 
lations shows that the easy axis of polarization is along 
[111] direction in BiFeOa. The ferroelectricity in BiFeOs 
is originating from the distortion of Bi-0 coordination 
environment as a result of the stereochemical activity of 
the lone pair on Bi. 

9. Large anisotropy in the polarization in BiFeOs has 
been identified and this could explain why there is large 
scattering in the experimentally reported polarization 

values in this material. 

10. Our partial polarization analysis demonstrates that 
the polarization is essentially due to the displacement of 
Bi relative to the center of the FeOg octahedra. 

The overall conclusion is that our results are in good 
agreement with experiment. This demonstrates that ac- 
curate density functional calculations can be used to ex- 



plore not only the structural and magnetic properties but 
also the ferroelectric properties of magnetoelectric mate- 
rials. 
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